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I. INTRODUCTION 



A model is proposed to explain temporal patterns of activity in a class of periodically exploding 

• Strombolian-type volcanos. These patterns include major events (explosions) which follow each other 
I every 10-30 minutes and subsequent tremor with a typical period of 1 second. This two-periodic 

^SJ ■ activity is thought to be caused by two distinct mechanisms of accumulation of the elastic energy 

in the moving magma column: compressibility of the magma in the lower conduit and viscoelastic 
2 , response of the almost solid magma plug on the top. A release of the elastic energy happens when 

a stick-slip dynamic phase transition in a boundary layer along the walls of the conduit occurs; this 
phase transition is driven by the shear stress accumulated in the boundary layer. The first-order 
character and intrinsic hysteresis of this phase transition explains the long periods of inactivity in 
the explosion cycle. Temporal characteristics of the model are found to be qualitatively similar to 
the acoustic and seismic signals recorded at Karymsky volcano in Kamchatka. 

^ ! 

^ ^. Keywords: stick-slip phase transition, viscoelasticity, hysteresis. 

o 

O ■ A wide variety of types of volcanic activity exists: from devastating explosions separated by calm periods of many 
' S ' years or even centuries to a steady continuous outpouring of magma. In this paper we study mechanisms that give 
rise to a Strombolian-type volcanic activity, which is somewhere between these two extremes and characterized by 
' regularly repeated explosions (10-400 per day) followed by relatively calm periods. Often in the course of a longer 
explosion, the gas and ash emission exhibits audible and visible modulations with a rather robust period of about 1 
second. These modulations, often found on the seismograms of Strombolian-type eruptions, received a special name, 
T— I [ "chugging" (Benoit and McNutt, 1997), because they resemble a periodic noise produced by a steam engine. An 
^ ■ example of a volcano exhibiting such activity, we consider Karymsky volcano on Kamchatka peninsula in the Far East 
' of Russia. 

We believe that these Strombolian-type eruption patterns arc caused by peculiarities of the motion of the magma 
column, and propose a model that describes the magma motion in a volcanic conduit as a creep flow of viscoelastic 
, compressible medium with shear-stress-dependent boundary conditions. We show that the dynamics of the model 

• indeed exhibits 2 levels of quasiperiodic behavior and resembles the eruption patterns of the Strombolian-type volcanos 
both on long and short time scales. The physical transparency of the suggested model allows us to express the observed 

' time scales using the material properties of the magma. 
O I Our paper is organized in the following way: first, we briefly describe geological aspects of the current eruption 
' 55 ' '^^ Karymsky volcano. Then, after referring to the existing explanations of Strombolian-type activity, the model is 
formally introduced. We derive simple analytical expressions for a dormant time, a duration of explosion , and a 
1 typical period of tremor, and then present a numerical solution to a system of dynamical equations describing the 
motion of the magma column. The paper is concluded with a discussion of results and possible directions for further 
studies. 

X 
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II. THE KARYMSKY ERUPTION OF 1996 - 2000 




Fig. 1. Karymsky volcano. 

Karymsky volcano (Fig. 1), located in the central part of the East Kamchatka volcanic belt, is one of the most 
active volcanoes in the far-eastern Russia. It is a typical andesite stratovolcano composed of a 7700 — 7800-year old 
caldera 5 km in diameter (Ivanov et al., 1991) and a more recent cone which started growing about 5300 years 
ago and at present (September 2000) is ~ 700 m high. The cone is composed of accumulated lava and pyroclastic 
materials, and the summit of the cone is capped by a double crater. The active part of the crater, formed during the 
most recent eruption, widened from 90 m in 1996 to 190 m in 2000 (Fig. 2). The absolute height of the volcano is 
1549 m. 




Fig. 2. Double summit crater at Karymsky volcano. 

Karymsky activity usually consists of long effusive-explosive eruptions, with the eruption previous to the one being 
studied taking place in 1970-1982. The current eruption at Karymsky volcano began in January 1996 and apparently 
entered a lull phase in December 2000. Each eruption consists mostly of the discrete quasiperiodic gas and ash bursts 
sending plumes 100-1000 m above the crater rim (Fig. 3). 
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Fig. 2. Gas and ash plums at Karymsky volcano. 

Ash and steam clouds have been observed extending downwind for 10-50 km, rarely 100-200 km. 

A lava flow field, produced by the very viscous magma, runs down the south-western slope of the cone and reaches 
the length of ~ 1.4 km and the width ^ 100 — 200 m. Even during the periods of highest activity, the motion of the 
lava flows was very slow (few meters per day) and was noticeable only by periodic rock-falls, happening every 5-15 
minutes. 

During particularly strong explosions, the eruptive column collapsed causing pyroclastic flows. In the lower part 
of the eruptive column, hot volcanic bombs of the sizes up to 2, rarely 5 m, were frequently observed. The bombs 
were dense and often covered with a bread-crust-Uke structure formed during cooling. On the surface of some bombs, 
noticeable extrusive tracks were found; they were scratched on the very viscous moving magma by the protruding 
parts of the magma channel. 

Total volumes of pyroclastic and effusive materials deposited on the cone of Karymsky during the eruption of 
1996-2000 is estimated to be 0.0274 and 0.0229 km^, respectively. Eruptive products, which include magma, volcanic 
bombs, and ash, predominantly consist of andesite, usually of black and dark grey color. The most common mineral in 
the eruptive products is plagioclase (20-25 %). Crystals of plagioclase can reach a size of 3-5 mm. The concentration 
of olivine and pyroxene is less than 1%, their crystal sizes are usually 2-3 mm. A chemical content of the andesite, 
averaged over 6 samples from the 1996 - 2000 eruption, is listed below: 5*102 - 61.54, Ti02 - 0.86, AI2O3 - 16.67, 
Fe203 ~ 2.43, FeO - 5.12, MnO - 0.11, MgO - 1.99, CaO - 5.31, TVaaO - 3.69, K2O - 1.58, P2O5 - 0.25. 

The feeding system of Karymsky volcano was studied by various methods, including gravity surveys, aeromagnetic 
surveys, photogrammetry, seismology and geodesy (Zubin, 1971; Maguskin, 1982; Shirokov, 1988). Despite shght 
discrepancies in the estimates of the depth of the upper border of magma chamber and its size, these data agree 
that underneath the Karymsky volcano in the close proximity to the surface there exists a spherical magma chamber. 
The upper border of this chamber is in few kilometers below the sea level. The diameter of the magma chamber is 
estimated to be between 1.5 and 7 km. A magma channel of a diameter of 100-200 m leads from the magma chamber 
to the top crater. 

III. DUAL PERIODICITY OF KARYMSKY ERUPTION 

A distinctive feature of Karymsky eruptions is its rather robust periodicity which is observed on two timescales. An 
explosion, and a following quiet period, define the first period of activity which varied between 3 and 20 minutes and 
sometimes extending up to an hour. An explosion itself can start either abruptly or gradually (Fig. 4) and proceed 
according to one of the following scenarios: 
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Karymsky Acoustic Boom-Chugs 




Fig. 4. Typical acoustic signal of a period of activity of Karymsky volcano 

during shorter explosions, the activity decays quickly and monotonously (in less than 10 seconds), while the longer 
explosions (20 — 40 s, rarely 1 minute) may contain periods of relatively high and low activity. Often during longer 
(> 20s) explosions, the intensity of the acoustic and seismic signals emitted by the volcano exhibits chugging, i.e. 
modulations with a typical period of ~ 1 second. The chugging defines the second, shorter, period in the cyclic activity 
of the volcano. A typical chugging event consists of 10-20 cycles (Fig. 4). The explosion is followed by a quiet period 
characterized by a complete lack of any activity in the crater and usually lasting much longer than the explosions 
themselves. Similar temporal patterns have been also observed during previous Karymsky eruptions. (Tokarev and 
Firstov, 1967; Farberov et al., 1983). 

The systematic seismic and acoustic studies of Karymsky volcano were conducted in 1997-1999 by 3 Russian-US 
expeditions. The seismograms presented in Fig. 4 were recorded in 1998 by a three-component wide-band seismometers 
CMG 40-T, but the most detailed information about a fine structure of the chugging signal was obtained using a set 
of infrasound microphones. For a more complete description of the seismic and infrasound recordings at Karymsky 
the reader is referred to (Johnson et al., 1998; Johnson, 2000; Johnson and Lees, 2000). 



IV. EXISTING MODELS OF ERUPTION PERIODICITY 



Currently there is no single well-established view on the mechanism causing periodic patterns in the eruptions 
of Strombolian type. Below we mention several possible explanations of origins of the periodicity. Lees and Bolton 
(1998) compared the eruption process to the steam exhaust from a hydrodynamic system with a nonlinear ly controlled 
exhaust valve resembling a pressure cooker. Manga (1996) suggested that the periodicity is caused by a segregation 
of ascending bubbles into waves, as in a pint of Guiness beer. Julian (1994) proposed a mechanism of nonlinear 
hydrodynamic oscillations in the ascending magma flow. Jaupart and VergnioUc (1989) suggested that oscillations are 
caused by periodical collapse of a bubble foam trapped in pockets of the magma chamber. Recently, there appeared a 
paper by Delinger and Hoblit (1999) where an ad hoc hysteresis model was developed to describe a periodic behavior 
of silicic volcanos. However, self-consistent explanation for both long- and short-time periodicity and extended calm 
periods in the activity of Strombolian-type volcanoes is still missing. Our model, based on viscoelastic hydrodynamical 
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description of the motion of the magma column with shear stress-dependent boundary conditions, accounts for both 
timescales. 



V. THE MODEL 

Let lis first qualitatively discuss the motion of magma in the conduit. Since the temperature and pressure in the 
magma column decrease with the height, the viscosity of magma increases rapidly in the upper part of the conduit. 
Inspecting eruption products ejected to the surface, it is natural to deduce that at the very top of the column the 
magma is almost solid. However, deeper parts of the magma column are still hot and much less viscous. Given this 
significant difference in rheological properties of the upper and lower sections of the magma column, we consider the 
moving magma as consisting of two distinct parts: long lower part (LP) and short upper part (UP). We assume that 
the viscosity of the magma in the LP is negligible. Due to the large amount of dissolved gases, magma in the LP is 
compressible. We also assume that there is a constant supply of fresh magma to the bottom of the LP. Unlike the 
LP, magma in the UP is cold and viscous; in addition to that, as any media near the liquid-solid transition point, 
it possesses a certain degree of elasticity. Because of the small length of the UP compared to the LP 100 m vs. 
~ 5 km), we can neglect compressibility of the UP. In summary, we consider a long cylindrical tube filled with a 
non- viscous compressible media, fed to the bottom with a constant velocity and a short viscoelastic plug (UP) on the 
top. The boundary conditions for the moving UP are controlled by a local value of a shear stress in the vicinity of 
the cylinder wall. When the boundary shear stress is low, the plug "sticks" to the cylinder wall and the boundary 
velocity is zero. As the shear stress increases, the magma near the wall undergoes phase transition (melts) to a much 
less viscous liquid state, and the boundary layer of the plug slips along the walls of the channel. This shear-induced 
phase transition happens when the viscoelastic state with the non-zero shear modulus becomes thermodynamically 
less stable than a pure liquid state with zero shear modulus. Phenomenologically, the transition between stick and 
slip boundary conditions for the very viscous body is similar to the effect of dry friction, where a motion begins only 
when a driving force exceeds some threshold value (Persson, 1999, Ch. 8). 

Qualitatively, a cycle of the system evolution can be described as follows: when the pressure in the compressible 
LP of the column is low, the shear stress (7 in the plug is much smaller than the threshold phase transition value g slip 5 
and the velocity of the UP is much less than the feeding velocity Vq- The difference in feeding and plug velocities 
results in the growth of the pressure of the compressible magma in LP which, in its turn, increases the shear stress in 
the plug. This stage corresponds to the dormant period between explosions. When the shear stress in the boundary 
layer of the UP exceeds the critical value asHp, the stick-slip phase transition occurs and the plug begins to move. 
This corresponds to the appearance of visible signs of a volcanic explosion. 

Depending on the parameters of the model, the further evolution of the system proceeds according to one of the 
following two scenarios. For certain conditions, a relaxation of the accumulated in the UP viscoelastic deformations 
can result in a low, or even reverse, velocity of the boundary layers of the UP relative to the walls of the conduit. If 
these conditions persist for sufficient time, they can cause a drop of the shear stress below the slip-stick transition 
critical value (TsUck, which in turn can give rise to the reverse, slip-stick transition. However, the stick state will be 
very short-lived, since the accumulated excessive pressure in the LP and oscillatory motion of UP will result in the 
immediate increase of the boundary shear stress and a new stick-slip transition. Length of a period of such coupled 
viscoelastic-stick-slip process is controlled by the density, the elastic modulus, and the characteristic time of the phase 
transition. In the second scenario, sticking does not occur during short-periodic viscoelastic oscillations of the UP, and 
the period of these oscillations is roughly equal to that of a freely oscillating elastic membrane with mixed boundary 
conditions. Either of these oscillation processes can cause the short-period modulation of the activity of the volcano. 

During the active phase, the velocity of the UP is much higher than the feeding velocity Vq which results in an 
expansion of the LP and a decrease of pressure under the UP. When the pressure under the UP drops so that the 
shear stress in the UP is below the slip-stick threshold value a stick, a "long-term" sticking of the plug occurs, and 
the system enters a new quiet period of evolution. Because of the hysteresis {(TsUp < GsUck) associated with the first 
order nature of the stick-slip transition, it takes some time to build enough pressure for a new explosion to begin. 
This cycle corresponds to the second, longer period (10-30 minutes) in the dynamics of eruption. 

Although simple in nature, this model qualitatively explains both scales of periodicity of Karymsky volcano and 
does not contradict to the observations. Below a formal analysis of the dynamics of the model is presented. 

Let us consider a cylinder (LP) of the radius R and the length L, filled with an ideal non-viscous compressible 
medium feedcd from the bottom with a constant velocity Vq. The top of LP ends with a short viscoelastic plug (UP) 
of length /, / < L (Fig. 5). 



5 



1 UP 




|v(r,t) 












LP 









Fig. 5. Sketch of the magma conduit. 



The velocity field in the UP, of density p and driven by pressure P and stress tensor a^j, is described by the 
Navier-Stokes equation, 



dP da, 



(1) 



Because of the high viscosity and the low velocity of magma in the UP, the nonlinear (convective) term is omitted. 
To account for the viscoelastic properties of the UP we use the simple linear Maxwell model (see, for example. Shore 
et al.), with a characteristic memory time r^. 



da. 



(2) 



Here 77 is the dynamic viscosity, and e^- 
the magma in the UP reads: 



^"^ dt ~ ' 
dvi/dxj +dvj/dxi is the shear rate tensor. Incompressibility condition for 



dvj 
dxi 



0. 



(3) 



For simplicity, we assume that the motion of the plug is axially symmetric and select the Z direction of our 
cylindrical coordinate system parallel the axis of the cylinder directed along the magma flow. Also, we also consider 
the pressure gradient constant through the UP and being equal to —j^, where Pq is a pressure at the bottom of the 
UP. Using (||J|), the Navier-Stokes equation (|^) for I4 (Z-component of the velocity of the plug averaged over the 
plug length) can be transformed to: 



dV,{r,t) Po{t) 



dt 



I 



— / cxp 

Tr 



t-t' 



A,K(r, t')dt'. 



(4) 



Here = /dr^ + 1/r d/dr is the radial part of the Laplacian operator in the cylindrical coordinates. 

From the conservation of mass of the compressible media in the LP, it follows for the pressure Pq on the bottom 
surface of the UP, 



v{t) - - y{t) ^ 2 



K(r, t)rdr/R\ 



(5) 



where (3 = —l/v{dv/dp)T is the compressibility (isothermal) of the media in the LP, and V{t) is the velocity of the 
UP Vz(r,t), averaged over the perpendicular crossection. 

To complete the description of the dynamics of the system, the equations (|^J^) must be complemented by boundary 
conditions for the velocity on the conduit walls, Vz{r — R,t). These boundary conditions are controlled by the kinetics 
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of the first-order phase transition driven by the value of the shear stress a{t) = arz{t, R) near the wall and are written 
in a general mixed form: 



V,{R,t) = -RiP{t) 



dr 



(6) 



r=R 



The quantity Rip{t) is usually called a "slipping length". It shows at what distance inside the wall the linearly 
extrapolated velocity becomes equal to zero. For = the boundary conditions for the velocity are "stick", i.e. 
Vz{R,t) = 0; for ip{t) > the boundary layers of the UP "slip" along the walls of the cylinder with some finite 
velocity, Vz{R,t) > 0. 

The time evolution of the dimensionless slipping length ip(t) is determined from the time-dependent Ginzburg- 
Landau equation, which describes relaxation of a system with a first-order phase transition: 



TGL 



dt 



dF 



(7) 



where tql sets a phase transition timescale. The free energy -F(V') is a usual double- well potential with a tilt depending 
on the value of the shear stress a{t) on the boundary of the UP. We choose a double-parabola potential which has 
fixed positions of both minima, i]j = Q and ip ^ ipQ: 



^0? 



tp < 1p' 
h, V > -0' 



(8) 



Here /i is a tilt parameter; matching point ip' is determined from continuity requirement, -0' = (V'o ~ h)/2ipQ. To 
determine h and to illustrate how Eqs. (^) work, let us consider a typical stick-slip cycle. We go around a hysteresis 
curve shown in Fig. 6. 




-*stick "slip 

Fig. 6. Hysteresis curve I-II-III-IV for the stick-slip and the slip-stick transitions. Inserts show typical sketches of F(%1)) for 
all branches of the curve. 



When the shear stress a{t) is small and growing (line I), the boundary condition is "stick" and stationary, -0 — 
^ = 0. As the shear stress exceeds the critical value CTsUp, ^{t) starts moving towards tjjQ with the rate ~ 2-00 /^gl 
(line II). After p reaches '00 j the boundary condition becomes "slip" and stationary again, and the shear stress starts 
decaying (line III). After a{t) drops below a stick, ip relaxes to (Hue IV). The state I corresponds to the dormant 
period, the state III corresponds to the explosion without intermediate sticking; the transition intervals II, IV are 
usually very short. The free energy tilt parameter h is linearly connected to a: 



2a 



^slip H" ^stick 



slip ^ stick ^slip ^ stick 



(9) 



Similarly to the Eq. (Q), in the framework of the Maxwell model the shear stress a{t) is expressed as a convolution 



of a shear rate 



dr 



with an exponential memory function exp(— -^^r-^) 



* , t-t\ dVz{r,t') 
exp( ^ 



dr 



dt' 



(10) 



r=R 
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The equations (^-[lO|) completely define the dynamics of our model. A very similar set of equation was derived by 
Shore et al. (1997) for the physically equivalent model of polymer extrusion. Because of the strong non-linearity of 
the Eq. an analytic solution of these equations is impossible. We refer readers to Shore et al., 1997, for a detailed 
description of a linear stability analysis of this system. Contrary to a stable flow, a dynamical regime corresponding 
to volcanic activity is strongly unstable and nonlinear. However, large separation of characteristic timescales allows 
us to give reasonable analytical estimates of important temporal features of the model. 

Let us first look at the shortest timescale which corresponds to small oscillations of the UP when it slips along the 
conduit walls without intermediate sticking. For this timescale the pressure term in is considered to be constant. 
Dynamics of the UP is physically equivalent to the oscillation of an elastic circle membrane with the mixed boundary 
conditions (p|). Only the lowest harmonic is considered, 



V,{r,t)=Voe-'^'Jo[^x) 



(U) 



where A, a number of order of one, depends on the slipping length ip: Jo (A) = Ji{X)X'ip. Here Jq and Ji are the zero- 
and first-order Bessel functions. After plugging Eq. (|ll|) into Eq. and disregarding the constant pressure term, we 
obtain the following dispersion relation: 



27 ■ 

This yields for a period Tqsc of not too overdamped UP oscillations {Tqsc ^ '^m)'- 



(12) 



Tos.^'-^^. (13) 

Taking into account a commonly used relation between the Maxwell time and shear modulus G, w 'y/G (see, for 
example, Webb (1997)), we observe that Eq. (^3|) reduces to a usual expressions for a period of a small free oscillations 
of an elastic membrane, 

2.R . . ^^^^ 



dV^ (r) 
dr 



Now we go to much longer timescales, T {Tosc,Tm} and look first at the steady state UP boundary shear stress 
. Since this timescale is much greater than the elastic memory time t^, the viscoelastic term in (0) 

r—R 

is relaxed to a usual viscous term, r/ArVzir,t). For a steady state, the averaged over the crossection of the conduit 
velocity V{t) must be equal to Vq. For an arbitrary slipping length ip, a straightforward calculation yields: 

_ Wpr, 1 

It immediately follows that in order for stick-slip repeating cycles to happen, the "stick" steady state shear stress {a 
for t/j = 0) should be greater than asUp and "slip" steady state shear stress (tr for -0 = V'o) should be less than a stick- 
These conditions imply for the feeding velocity Vq: 

<ysiipR o-sticfe-R(l + 4i/'o) 

< ^0 < ^ . (16) 

Besides the steady state properties, we can also evaluate characteristic times between the stick-slip and the slip- 
stick phase transitions, i.e. duration of the dormant and explosive periods. Assuming that these times are much 
longer than both T^sc and t^, we again disregard all memory effects in the viscous term in Eq. (^ and consider the 
flow inertialess by omitting the left-hand-side term After averaging the remaining part of Eq. (^ over the 

crossection of the flow and plugging it in (^, we obtain for arbitrary slipping length ■0: 



/' 



I3LI i?2(l-t-4V') 

Combining the solution of this equation, 



''"-^'('\,_3M2L.o. (17) 
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Vit)^K + Ce-^, to^j^^^y (18) 



with expressions for stick and slip velocity that follow from (|15D, we obtain for dormant (Td) and explosive (Tg) times: 

^ _ 8tj(3LI Af]Vo - asUckR , . 

' R^il + i^jp) ^^asUckRil + ^M-^vVo' ^ ' 



VI. NUMERICAL RESULTS 

To check our analytical predictions and get an overall view of the system dynamics, we solved the Eqs. (|^-[To|) 
numerically. We used a simple Eulerian finite-difference scheme with a uniform radial grid of 50-100 points. All the 
integrals are approximated by trapezoid formula. Given that time step is sufficiently small to avoid von Neumann-type 
instabilities, this simplest difference scheme proved to be sufficiently robust. 

We chose the numerical parameters to be in general accordance with the literature (Murase and McBirney, 1973;, 
Borgia and Linneman, 1990; Griffiths and Fink, 1993;, Bagdassarov and Dorfman, 1994; Webb, 1997; Bagdassarov 
and Dorfman, 1998) and with our estimates ( |l^JT^ , pj| , pO| ) . The following values were used: 

Density p = 2400 kg/m^ 
Dynamic viscosity ry = 10® Pa s 
Conduit raidus R — 50 m 
UP length I ^ 100 m 
LP length L = 5 km 
Compressibility /3 = 4 10"^° Pa"! 
Maxwell time r™ — OAs 
Ginzburg-Landau time tgl — 10^'^ s 
Average velocity Vq = 2 x 10^^ m/s 
Slipping shear stress asUp — 10^ Pa 
Sticking shear stress a stick = 310"* Pa 
Dimensionless slipping length ?/; = 10 



With this choice of numerical parameters a no-stick explosion takes place, during which the UP moves along the 
walls of the conduit with the constant slipping length 'i\)R until it finally sticks to the wall at the end of the explosion. 
For such explosions a choice of phase transition time tgl does not affect the dynamics, given that it is shorter than 
explosion time, tgl ^ Te. 

Plots of averaged over crossection velocity V(t) and pressure gradient vs. time are presented in Fig. 7. The 
simulation results for dormant and explosion times, as well for short-term oscillation period, « 530 s, « 21s, 
T^osc ~ 0-69 s are in good agreement with our simple analytical estimates, Td — 497 s, Te — 20 s, and Tosc — 0.69 s. 
To calculate Tosc we have to solve the equation Jo(A) = Ji(A)A'0, which for -0 = 10 yields A « 0.43. For the stick 
boundary conditions {ip = 0), A w 2.3; which gives for the period of oscillations at the end of explosion cycle when 
UP sticks to the walls T^^^, « 0.13, which is again in very good agreement with the simulation data. 
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Fig. 7. Time dependence of the spatially averaged UP velocity (solid line) and the rescaled pressure gradient (in 5 x 10'' Pa/m, 
dashed line). 

Another possible scenario for an eruption is when during explosion cycle the UP sticks to the walls numerous times. 
It happenes when the shear stress at the wall, modulated by the free membrane-like oscillations of the UP, drops 
below the sticking value dstick, and the slip-stick phase transition is sufficiently fast, tgl *C Tqsc- To simulate this 
regime, we increase the value of aguck from 3 10'* to 8.5 10^ Pa and keep all other parameters constant. The long-time 
dynamics, characterised by the times Tg and only slightly changes since both of these times depend on agUck 
logarithmically. However, a short-time dynamics, previously described by decaying free oscillations with a period 
Tosc, changes radically and becomes chaotic (Fig. 8). 
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A typical time separation between consecutive velocity maxima (w 0.3 s for the first three maxima in the bottom 
part of Fig. 8) becomes somewhat smaller than free oscillation period Tosc « 0.7. Qualitatively, it happenes when 
during an oscillation the velocity increases so that shear stress drops below u stick, a sticking occurs that results in a 
rapid decrease of velocity, so the oscillation is "truncated" before it reaches its maximum. Another distinct feature 
of this regime is that the short-time chaotic modulation of the UP velocity continues through the whole explosion 
phase, while in the non-sticking regime, the oscilations decay after a few periods. 



VII. DISCUSSION 



In the previous sections the long-term behavior and two possible scenarios of the short-time activity of the model 

based on viscoelastic hydrodynamics and stick-slip transition were described. Patterns in the time dependence of 
the UP velocity, presented in Figs. 7,8, qualitatively match the seismic and acoustic signals (Fig. 4) recorded at 
Karymsky volcano. Yet nothing has been said about a mechanism describing generation of seismic and acoustic waves 
themselves. The oscillating UP and related pressure oscillations in the LP can cause seismic waves either directly or 
act as a "pacemaker" for other processes. These processes can include a stimulated emission of gases either through 
creation and change of geometry of cracks in the UP or via degassing of magma under action of standing and running 
compression waves. Since these mechanisms have a strongly nonlinear character, their coupling to the UP oscillation 
can modify the oscillation period and decay patterns. We leave consideration of these phenomena for future study. 
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VIII. CONCLUSION 



As a result of the visual, seismic, and acoustic observations at Karymsky volcano, two scales of temporal periodicity 
were observed in the eruptive dynamics. First, the repetitive explosions which produced tha ash and gas clouds and 
volcanic bombs in 10-20 minutes intervals, and secondly, the periodic modulation of eruptive activity (chugging) with 
a typical timescale of order of 1 s. The existence of these two scales of temporal periodicity are not limited to the 
1996-2000 eruption of Karymsky volcano, similar timescales were observed during previous Karymsky eruptions and 
on other andesite volcanos in the world. 

A model of the motion of the magma column was suggested that accounts for both timescales in the eruption of 
Karymsky and other similar andesite Strombolian-type volcanos. Low- and high-frequency periodicity is caused by 
two distinct ways of accumulation of the elastic energy. The viscoelastic properties of the UP causes the high frequency 
tremor, while the compressibility of magma in the LP explains the low frequency cycles. The shear-stress-dcpendcnt 
stick-slip phase transition in the UP introduces the hysteresis into the dynamics of the magma column motion and 
explains the unharmonicity of the oscillations and long dormant periods in the eruption dynamics Physical simplicity 
of the model allows analytic estimates for the long and short periods of the explosive activity to be obtained. 
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